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, By discussing the Cauchy problem, we determine the covariant equation of the 

characteristic hypersurfaces in a relativistic superfluid theory. 
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^ 1 Introduction 

(2J[). Phenomenological approaches to relativistic superfluidity had been initially proposed by 
Rothen [Q, Dixon and Israel || ||. Then, Khalatnikov and Lebedev || and Carter 
|7], H have developed general superfluid theories which are in fact equivalent as they have 
later shown || |K| . This is a relativistic generalization of the non-dissipative version of the 
Landau theory of superfluidity. 

In the framework of this superfluid theory, Carter and Langlois [ITJ have studied the 
velocity of propagation of the sound following a method due to Hadamard. The story 
of infinitesimal disturbances, giving the sound speed, is described by an hypersurface for 
which the field variables are continuous but their derivatives are discontinuous across this 
one. The problem is to find the equation of these characteristic hypersurfaces. However, 
Carter and Langlois do not solve the general case but the low-temperature limit for 
a relativistic superfluid with phonon-like excitation spectrum whose they have derived the 
Lagrangian. Likewise, Vlasov |12|, [13] analyses the shock waves in this relativistic superfluid 
theory. 

The purpose of the present work is to give the covariant equation determining the 
characteristic hypersurfaces in this relativistic superfluid theory. Our method is based 
on the fact that the Cauchy problem of the field equations has no unique solution on a 
characteristic hypersurface. To do this, we consider the covariant field equations of the 
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superfluid in a spacetime with a background metric g pu in a general coordinate system 
We attempt to solve the initial value problem on the hypersurface x 1 = by finding 
the Cauchy data and by expressing merely the derivatives with respect to x 1 in terms of 



the Cauchy data, method used in hydrodynamics for instance by Lichnerowicz | l3fl . If it 
is not possible to determine uniquely the derivatives with respect to x 1 then x 1 = is 
a characteristic hypersurface. Also, discontinuities of the normal derivatives of the field 
variables can occur across x l = 0. In another coordinate system (y 11 ), the equation x l = 
will be f(y p ) = and by this way we will find a covariant equation for the function /. 

The plan of the work is as follows. In section 2, we recall the field equations of the rela- 
tivistic superfluid theory. We discuss the Cauchy problem in section 3. The determination 
of the covariant equation of the characteristic hypersurfaces is carried out in section 4. We 
add in section 5 some concluding remarks. 



2 Field equations of the relativistic superfluid 

The relativistic superfluid theory of Carter is built on a particle number n p and an entropy 
current s p which are conserved. In a general coordinate system, these conservation laws 
can be written 

vx = o , (i) 

V p s p = (2) 

where V CT denotes the covariante derivative associated with the metric g pv . A master 
function A of n p and s p defines the chemical potentiel \i a and the temperature 9 a by the 
differential relation dA = fi p dn p + 9 p ds p . The convective variational principle of Carter 
yields the equation 

s a (vj p - v p e a ) = o (3) 

and a similar equation for \i a but the assumption of superfluidity requires the compatible 
condition 

V CT/ u p - V p /v = . (4) 

For our purpose, it is more convenient to introduce a Lagrangian L having the expression 
L = A — fi p n p which is considered as a function of s p and \i a defining Q a and n p by the 
differential relation dL = 9 p ds p — n p dfi p . The Lagrangian L is in fact a function of the 
following scalars 

s 2 = —s p s p , /i 2 = —fi p fi p and y 2 = —s p fi p (5) 

defined for a signature ( — h ++) of the metric. In consequence, we can express n p and 9 P 
in terms of \i a and s p in the following form 

n p = Bfi p -As p and 9 p = Cs p + Api p (6) 
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where the functions A, B and C of s 2 , /i 2 and y 2 are the partial derivatives of L 



» = . C = -2^ and A=-£. (7) 



We henceforth consider that the field variables of the relativistic superffuid theory are 
jjjfj and s p . From ([!]) and (|3|) we derive the following equations 

/2 p d p B + BWy - s p d p A = , (8) 

sW , G + Cs^ + + + i<W + = (9) 

where we have used (|2]) and The derivatives of A, B and C appearing in (|) and (|) 
have the expressions 

9a A = ^d a s 2 + |^ 2 + |^y 2 . (10) 

We notice the identities 

dA IdB dA _ 1 <9C <95 dC 

d/j, 2 2 ch/ 2 ' 9s 2 2 <9y 2 9s 2 <9p 2 

For a Lagrangian L function of s 2 , /i 2 and y 2 , the field equations are (0), (f|), (§) and 
(|) with definitions (0) and (0). 



3 Cauchy problem 

We now discuss the Cauchy problem on the hypersurface x l = for the field equations of 
the previous section, the field variables being \i a and s p . The Cauchy data are s p and \i a 
on x 1 = 0. The problem is to determine d\\i a and d\s p on x 1 = in terms of s p and \i G 
on i 1 = and their derivatives with respect to x a (a = 0,2,3). Hereafter, we call dC all 
quantities calculable in terms of the Cauchy data. 

We immediately remark by virtue respectively of (0) and (^) that 

d lfJla = dC (a = 0, 2, 3) and d^ 1 = dC . (11) 

According to (§), we thus have 

«9 ly u 2 = -2 / u 1 9 1 /i 1 + dC and d lS 2 = -2s a d lS a + dC . (12) 

We insert ([TT|) and ([TJ) in equations (§) and (|9|) for p = a and we obtain four linear 
equations determining <9i/ii and <9is a in terms of the Cauchy data 
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The hypersurface x l = is characteristic if the determinant of the matrix in ([13]) van- 
ishes. Unfortunately, it is not easy to find the covariant equation of the characteristic 
hypersurfaces in this manner. 

Instead of making use of (0), we are going to show that equations @ and (H) determine 
difii, dis 2 and d\y 2 in terms of the Cauchy data 

/ mu ™>12 m 13 
777.21 ^22 "m 23 

\ 77731 m 32 m 33 
Furthermore, we will get 





s l d lS a = dC 

when difii, dis 2 and <9i?/ 2 are known from equations (|i4|). 
We have always equation (§) which gives 

l^dtB + Bg^dm! - s 1 d x A = dC . 

We add equation @ for p = 1 in the form 

(s 1 ) 2 ^ + sV^iA + s 2 g 11 d 1 C + ^Cg u d lS 2 + y 2 g 11 d 1 A = dC . 
taking into account ([□]). By contracting equation (§J) by fi a , we obtain thereby 



(14) 
(15) 

(16) 
(17) 



-y 2 s' 7 d a C-Cs a d a y 2 -Cs a s p V a n P -fi 2 s a d a A + s 2 n P d p C+-Cfi p d p s 2 + y 2 fi p d p A = (18) 

where we have used the identity 

s'UpVvS" = s°d a y 2 - s°s p V ai i p . 
We now write down equation (|T^) in the form 

-y 2 s 1 d 1 C-Cs 1 d 1 y 2 -C(s 1 ) 2 d 1 fi 1 ~fi 2 s 1 d 1 A+s 2 fi 1 d 1 C+^Cfi 1 d 1 s 2 +y 2 fi 1 d 1 A = dC . (19) 

In consequence, the three equations flTB"]), ( |i~?D and ( |T9"D give a system for determining 
d\H\, dis 2 and d\y 2 in terms of the Cauchy data . We obtain explicitly the following system 



-2(/i 



1\2 



OB 
dfi 2 



2 S y 



dA 
dfi 2 



+ Bg 



li 



<9i/ii 



+ 



x 55 x dA 

ds 2 ds 2 



d lS 2 + 



' 1 dB _ x dA 
^ dy 2 dy 2 



d lV 2 = dC , 



(20) 



' 1 2 d ° 2(/) - 2»V^s 2 - 2^V^y 2 



9/i 2 



+ 



,1\2 



dC 



a..i9A 



dfi 2 
ndA, 



<9i/i! 



<9C 

(* 1 ) 2 f§wl4 + 9 11 & 2 + 9 11 f^' 2 



1 



Cg 



n 



dy 2 



dy 2 



dy 2 



2 

<% 2 = dC , 



(21) 
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dfi 2 



+ 



,sc 2 ,dA 2 ,ac 2 , 1 i 

- S ^ + »d? S + 2 C » 

1 dC 2 iM 2 2 1 dA 2 „ J 

a?T a?/' 2 ajr a?r 



y 2 - C(« 1 
«9 lS 2 

9iy 2 = c/C 



(22) 



which corresponds to system (fH]). Moreover, from equation (|S|) for p = a we find ( P~5|) 
when the equations (EDI), (0) and (p2|) are satisfied. 



4 Covariant equation of the characteristic hypersur- 
faces 

When the condition detm^ = is verified it is not possible to determine uniquely <9i/zi, 
dis 1 and d\y 2 in terms of the Cauchy data. This condition defines the characteristic 
hypersurfaces. 

In a coordinate system (?/*) , for instance the Minkowskian coordinates, the equation of 
an hypersurface is f{y^) = 0. If we perform the change of coordinates 

x 1 = f( y ») and x a = y a (a = 0, 2, 3) (23) 

then the equation /(?/') = becomes x l = 0. We can then apply the results of the previous 
section in the coordinates (x^). It is crucial to notice that 

s 1 = s a d a f , f t L = f j, a d a f and g ll = g aP d a fdpf (24) 

where here s a , /i a and g a @ are the components in the coordinates (y^). Furthermore, s 2 , 
p 2 and y 2 are scalars which can be written in the coordinates (y M ). 

Now matrix m^- depends only on s 1 , /i 1 , g 11 and the scalars s 2 , fi 2 and y 2 taking into 
account the expressions of the functions A, B and C, therefore the condition det m^- = 
gives a covariant equation for the function / 

f(s a d a f, fi a d a f, g^djdpf, s 2 , /i 2 , y 2 ) = (25) 

which is the desired equation of the characteristic hypersurfaces. According to ([151) we 
must add the another possibility 

s a d a f = . (26) 
The speed of the sound u with respect to the superfluid fi p is given by 

(/^W) 2 

(jj?gP° + lififi^drfdrf 

in units for which c = 1. 
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5 Conclusion 



Of course, it is difficult to exploit equation (|25|) determining the characteristic hypersur- 
faces. For obvious reasons of causality, we must have u < 1 and so these characteristic 
hypersurfaces must be timelike, i.e. 



g aP d a fd p f > o . 



(27) 



By solving (p5|), we must find condition ( f2~7j ) for all solutions to the field equations corre- 
sponding to a given Lagrangian L. This restricts the generality of the function L. 

We can apply our formalism in the case of the Lagrangian found by Carter and Langlois 



11] which has the expression 



3k 
T 



2/3 



(28) 



where k is a physical constant and cp the velocity of sound for the fluid with one constituent 
defined by L such that cp < 1. We assume that the relative velocity v between \i p and s p 
defined by 

,2 



1 -v z 



2/ 4 

is smaller than cp. The temperature T is given by 



y 



.2 



1 2/3 



With expression (|2"8"|), we can evaluate the functions A, B and C and their derivatives that 
we express as functions of T, /z and v. We take the low-temperature limit, T — > 0, keeping 
/i and f fixed. Hence the condition det = reduces to two separate conditions 



lim mii = , 

lim T(m 22 m 33 - m 23 m 32 ) 







(29) 
(30) 



Condition (|29| ) gives immediately uj = cp. Condition (|30| ) is complicated and we are not 
enabled to obtain a simple formula in order to compare with the one of Carter and Langlois 
nj . For the case in which v = 0, we find un = cp/\/3 as expected. 
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